Diffusion-limited deposition of dipolar particles 
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Deposits of dipolar particles are investigated by means of extensive Monte Carlo simulations. We 
found that the effect of the interactions is described by an initial, non-universal, scaling regime 
characterized by orientationally ordered deposits. In the dipolar regime, the order and geometry 
of the clusters depend on the strength of the interactions and the magnetic properties are tunable 
by controlling the growth conditions. At later stages, the growth is dominated by thermal effects 
and the diffusion-limited universal regime obtains, at finite temperatures. At low temperatures 
the crossover size increases exponentially as T decreases and at T = only the dipolar regime is 
observed. 
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The growth of deposits by irreversible aggregation of 
particles is of great technological importance as well as of 
theoretical interest. A variety of mechanisms are involved 
in the growth processes but at late times scaling laws de- 
pending only on a few parameters have been observed. 
A general assumption, that describes the patterns found 
in many experiments, is that the deposition process is 
dominated by thermal diffusion. A simple model for this 
type of growth is diffusion-limited deposition (DLD) [l|, 
characterized by the formation of branched, fractal struc- 
tures similar to those found in electrodeposition, dielec- 
tric breakdown, etc (see ,21 and references therein). Un- 
der certain circumstances, however, interparticle inter- 
actions, favoring ordered structures that compete with 
the randomness of the diffusion process, are required to 
describe the observed growth patterns. A case in point 
is the diffusion-limited deposition of magnetic particles 
subject to dipolar interactions. 

Dipolar interactions are essential in determining the 
rich variety of structures exhibited by magnetic materi- 
als |3|, and their interplay with thermal diffusion may 
lead to novel magnetic properties. On the theoretical 
side, dipolar interactions provide a simple model to study 
the effects of anisotropic, long-ranged interactions on 
far from equilibrium aggregation processes. The cen- 
tral question concerns the change in the fractal dimen- 
sion of the aggregates, D, as the dipolar interactions are 
switched on. Meakin et al. 4] considered the effect 
of isotropic long-ranged, l/r*^, interactions in reaction- 
limited cluster-cluster aggregation (CCA) models. They 
found that D is unchanged for short-ranged interactions, 
i.e. for e > 2Do where Do is the fractal dimension of the 
non-interacting aggregates, while for longer-ranged in- 



teractions D may change substantially. Accordingly, nu- 
merical results for diffusion-limitted aggregation (DLA), 
performed for particles with Ising spins (short-range in- 
teractions), revealed no changes in the fractal dimension 
of the aggregates with increasing exchange interactions 
[5j. Finally, results for DLA 6] of dipolar particles in- 
dicate that D decreases as the strength of the dipolar 
interactions increases, in line with the results for CCA 
of dipolar particles 13| and with experimental results for 
the aggregation of magnetic microspheres ;;8| , but in dis- 
agreement with preliminary results of ours (9J . 

In this article we report results of extensive Monte 
Carlo simulations that provide a general framework 
where the apparently contradictory results described 
above may be understood. We show that the initial 
stage of two dimensional dipolar DLD growth is indeed 
described by a new nonuniversal scaling regime, charac- 
terized by clusters (trees) whose shape and fractal di- 
mension are temperature-dependent. For large enough 
systems, however, the dipolar regime crosses-over to the 
diffusion driven universal regime, where the effect of the 
dipolar interactions is dominated by thermal effects. It 
is also shown that the dipolar regime corresponds to ori- 
entationally ordered deposits and that the onset of the 
universal regime coincides with the disappearance of the 
orientational order. At T = only the dipolar regime is 
observed. 

In the new dipolar regime, the orientational order as 
well as the shape and fractal dimension of the clusters 
depend on the strength of the interactions. Thus, the 
magnetic properties of dipolar deposits may be tuned by 
controlling the growth conditions, such as temperature. 
Finally, we found that the fractal dimension of the entire 
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FIG. 1; (Color online) Average height H{s) and width W{s) 
as a function of the size of the tree, s, at T* = 10~^ and 
various system sizes. The lines are power laws |(2j with the 
exponents listed in table I. Dashed line: universal regime. 



FIG. 2: (Color online) Average height H{s) and width W{s) 
as a function of the size of the tree, s, at T* = 10~^ and 
various system sizes. The lines are power laws with the 
exponents listed in table I. Full line: dipolar regime. Dashed 
line: universal regime. 



deposit is always given by the universal (diffusion driven) 
value implying that in the dipolar regime, the trees have 
a fractal dimension Dt that differs from that of the entire 
deposit. 

We consider a two-dimensional box of side L and height 
H, on a square grid of side a = 1. Periodic boundary 
conditions are applied parallel to the substrate, of size L, 
and particles are dropped from a finite height above it. 
The particles carry a dipole moment of strength /z and 
interact through the pair potential 
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3(Ai •fi2)(A2 •i"i2) - Ai • A2 



(1) 



where ri2 is the distance between particles 1 and 2, fi2 
is the two-dimensional unit vector along the interparti- 
cle axis, and Ai ^-nd A2 ^^^ the three-dimensional unit 
vectors in the direction of the dipole moments of par- 
ticles 1 and 2 respectively. A particle is released at a 
height Hin with a dipole moment oriented at random. 
The particle undergoes a random walk through a series 
of jumps to nearest-neighbor sites, while interacting with 
the particles in the deposit. At each step a new po- 
sition and a new random 3-dimensional dipole orienta- 
tion are accepted according to a simple Metropolis rule 
based on the difference between the dipolar energies of 
the two configurations and defining the effective temper- 



ature T* — kBTa^/fi^. The long range of the dipolar 
interactions is taken into account by an Ewald summa- 
tion for the slab geometry of the system |9( . In the limit 
T* ^ 00 all displacements are accepted and the model 
reduces to DLD. A particle will eventually (i) contact 
the deposit or the substrate sticking to it irreversibly as 
its dipole relaxes along the local field [lOJ; or (ii) reach 
a height greater than ifout, when it is removed and a 
new particle released. See [9|, [llj for the details of the 
simulation. 

Simulations were carried out at 4 temperatures, 
T* = 10-\ 10-2, 10-3, 10-4, and 4 system-sizes, L = 
200, 400, 800, 1600 with 20, 30, 50, and 100 thousand par- 
ticles per deposit, respectively. The deposits are similar 
to those of DLD: they consist of many small trees com- 
peting to grow |9|. As the number of particles in the 
deposit increases, fewer and fewer trees keep on growing 
due to shadowing until only a single tree survives. 

We start by investigating the dependence of the height, 
H, and the width, W, of a tree with its size s (number 
of particles), as well as the distribution of trees Ug, i.e. 
the average number of trees of a given size. In DLD the 
trees scale as 0, llS ll^ , 
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FIG. 3: (Color online) Average number of trees of size s di- 
vided by L {fis/L) as a function of s, at various temperatures 
and system sizes. Lines as in figure 2. 



for sufficiently large values of s. Owing to the finite size 
of the simulation box, however, only a single tree survives 
when the number of deposited particles is large enough 
and, thus at large s the width of the tree saturates {W « 
L), then the height grows linearly with its size, {H oc s), 
and Hs exhibits a discontinuity. 

The difference between vu and p± measures the 
anisotropy of the trees. If vh = iy± the trees are isotropic 



and their fractal dimension is Dt = 
anisotropy is present, Dt becomes [^ 

A = 1 + (1-/^11 )/i^±. 



l/i^l|, whereas if 



(3) 



We note that the assumption that the deposit has the 
same fractal dimension as the trees, D = Dt, which holds 
in DLD, is not warranted in general. D was estimated 
through the average particle density at height h, p{h), 
which was found to scale as in DLD: at early times the 
deposit builds up until it reaches a height hi. Then, we 
found a scaling regime where the density decreases as a 
power law of the height, p{h) ~ h~°', with D = 2 — a [2|. 
The density saturates when the lateral correlation length 
^11 reaches the size of the system, at hg ^ U' . Simple 
arguments show 2] that the exponents vn , t and a are 
related through 
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(2-r)/^/||. 
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We have verified (jJl by calculating the average maxi- 
mal height and width of the largest tree, H{s) = {h{s)), 
W{s) = {w{s)), as a function of the tree size, s, and 
fis = (ns)- h{s) and w{s) are the maximal height and 
width of the largest tree of each deposit, and (...) is an 
average over all the deposits. 

Figures 1 and 2 show the results for H{s) and W{s) 
obtained at T* = 10"^ and T* = lO'^ for different box 
sizes. The points are histograms on a logarithmic scale. 
As expected, the results at T* = 10^^ correspond to 
those of DLD. Scaling was found for a wide range of s 
(s > 10), after a crossover from the early stage regime 
(s < 10). The crossover to the linear regime (only-one- 
tree growth) was observed only for the smallest system, 
L — 200, for s > 5000 (not shown). nw|| and vx calcu- 
lated for each L were found independent of system size; 
a representative value was obtained by averaging over all 
L with the uncertainty estimated as the largest devia- 
tion from the mean. Power law regressions for different 
ranges of s, yield i^n — 0.64(1). The crossover to satura- 
tion of the width was observed for sizes up to i = 800. 
Fitting only points that exhibit clear data collapse yields 
i'± — 0.60(2). These exponents agree with previous re- 
sults for DLD [l2|, and show that at high temperatures, 
the universal behavior of the geometrical properties of 
the trees is unaffected by the dipolar interactions. 

At T* = 10""^, however, the existence of two scaling 
regimes is apparent, in particular for the largest system, 
L = 1600, where the width saturation and the linear 
regimes are not observed. For s less than the crossover 
size s* « 500 we found i^|| = 0.74(1) and iy±_ = 0.78(2). 
These exponents differ from those of DLD and character- 
ize a new growth regime that we call the dipolar regime. 
For s > s*, we obtained, iy\\ = 0.64(1) and v±_ = 0.60(2), 
in line with the results for DLD. Similar behavior was ob- 
served at the other two temperatures, that is, a dipolar 
regime with temperature-dependent exponents followed 
by a second scaling regime with DLD exponents. 

In table I we list the exponents of the dipolar and uni- 
versal regimes, obtained at various temperatures. These 
results indicate that, if ^ is assumed, the effect of the 
dipolar interactions may be described by the appearance 
of a dipolar regime characterized by non-universal expo- 
nents i^ii and i/_L that increase with decreasing temper- 
ature. The crossover to the universal regime occurs at 
tree sizes that increase as the temperature decreases. As 
a consequence, the universal regime is difficult to observe 
at very low temperatures. However, at finite T*, the uni- 
versal regime may be reached if large enough deposits are 
grown. Indeed, we found no dependence of the crossover 
between the dipolar and the universal regimes on L, and 
thus at any temperature there is an L above which this 
crossover may be observed. This implies that the geomet- 
rical properties of the trees may be tuned by controlling 
the dipolar interactions and the system size. In particu- 
lar, it is possible to deposit trees with a given anisotropy 



by growing trees (at a fixed temperature) with a fractal 
dimension that is determined by the tree size. Alterna- 
tively, at fixed tree size, one may control the anisotropy 
by changing the temperature. 

In figure 3, we plot the results for the tree distribu- 
tion at T* = 10~^ and T* = IQ-^. At T* = 10-\ the 
results are those of DLD. Scaling behavior was observed 
for a range of s between 10 and the maximum tree size 
(which depends on L), with an exponent r = 1.54(2) in 
line with results for DLD ^^. At T* = 10"^, assuming 
the crossover to occur at s* « 500, we estimated r in the 
ranges 10 < s < s* (dipolar regime) and s* < s < 5000 
(universal regime). We found 1.40(2) and 1.56(4), re- 
spectively, for both systems. Thus, as for the tree height 
and width, we found a dipolar regime for the tree dis- 
tribution, characterized by a temperature-dependent ex- 
ponent. Again, the crossover to the universal regime oc- 
curs for sizes that increase as the temperature decreases. 
This crossover was not observed at T* = 10^^: due to 
the limited number of deposits and decreasing density, 
the number of trees larger than s « 500 — 1000 was too 
small to be analyzed. Thus, at this temperature only the 
dipolar regime was observed, with r = 1.36(3)). 

We have used © and Q to verify the consistency 
of the exponents and to estimate Df. In addition, we 
have checked the validity of Q) in the dipolar regime and 
compared Dt in both regimes, a was estimated in both 
regimes by power law regressions of p{h) in the regions 
suggested by the plots of Hi (using data from L = 1600). 
We found weak crossovers at the temperatures T* = 10~^ 
and 10^'^, similarly to what was observed for n^. The 
results listed in Table I indicate a remarkable consistency 
between the values of r obtained from simulation and 
using Q). It is also clear that the fractal dimension of 
the trees in the dipolar regime decreases with decreasing 
temperature, in line with previous studies of DLA dipolar 
aggregates laLDl- 

The connection between the orientational order of the 
dipoles and the geometrical properties of the deposits 
was investigated by comparing the mean-square magne- 
tization density in the x— and y— directions at height 
h, (rn^ y{h)), with the mean particle density at the same 
height, p{h). These results are shown in figure 4. It is 
apparent that the density saturates at the same height 
as the mean-square magnetizations, and that the onset 
of the scaling behavior of p coincides with that of (my). 
Notice also the height /i^ w 300, where (m^) — {'m't) 
and the orientational order vanishes. That h) signals 
the disappearance of orientational order is seen most 
clearly in the inset of figure 3. The latter shows that 
the diagonal elements, Qxx{h) and Qyy(h), of the order- 
ing matrix Q become identical at h = h'' , beyond which 
Qxx — Qyy ~ 1/4 [ll|. We also found that h^ does not 
depend on L, but increases with decreasing T* . At low 
T*, the increase is exponential, /i^ ~ exp(l/r*). This 




FIG. 4: (Color online) Average particle density and mean- 
square magnetization densities at height h. Inset: diagonal 
elements, Qxx and Qyy, of the ordering matrix Q. 



functional dependence was obtained from a direct esti- 
mation of the heights at which the orientational order 
of the deposits vanishes. Depending on the tempera- 
ture and system size, h^ is larger or smaller than hg and 
this is related to the existence of the dipolar and univer- 
sal regimes: when h^ < hg both the universal and the 
dipolar regimes occur, but only the latter is found when 
Ji^ > hg. We recalculated the scaling exponents assum- 
ing that the crossover between the dipolar and universal 
regimes occurs sA, h = h^ (equivalently, at s^ estimated 
using H{s)) and found that they are identical with those 
of table I. Therefore, the dipolar regime is characterized 
by deposits with orientational order. As a consequence, 
and by contrast to what happens in DLD, the fractal 
dimension of the deposits depends both on the fractal 
dimension of the trees, Dt, and on the distribution of the 
trees, r. In the light of this, the decreasing Dt{T*) re- 
ported in |3,|8l, is to be expected, since it was determined 
from the radius of gyration versus the (small) number of 
particles of each cluster. Furthermore, the results of p, 
obtained for a single DLA cluster, may be interpreted as 
the crossover from the temperature-dependent fractal di- 
mension at short length scales, to 13 ~ 1.7 at long length 
scales p. 

We have also estimated the average interaction of a 
dipolar particle with similar particles in ordered and ran- 
domly oriented deposits numerically, and confirmed that 
for ordered deposits the interaction decays more slowly 
than 2Do(= 2 x Dt{T = oo) k. 3.28) while for random 
ones it decays faster, in line with the results of |j| for 
isotropic systems [ijj. Finally, we have checked that the 
global fractal dimension of the deposits is unaffected (or 
very weakly affected) by the dipolar interactions. This 
was found by analyzing the decay of the mean density 
with h, the scaling of hg and of the density at saturation 
with L, the two-point density-density correlation func- 
tion, the initial divergence of the interface width, and 



TABLE I: Characteristic exponents of the dipolar DLD model, 
obtained using the equations indicated in brackets. 



and T obtained from the simulation; other exponents 
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the mean height of the upper surface. In every case no 
significant deviation from DLD was found [lj|. 

Funding from the Portuguese FCT, contract 
SFRH/BPD/5654/2001, and from the Spanish MCyT, 
project BMF2001-2841, is gratefully acknowledged. 



[1] 
[2] 



[3] 
[4] 
[5] 
[6] 



P. Mealiin, Phys. Rev. A 27, 2616 (1983). 
P. Meakin, Fractals, scaling and growth far from equilib- 
rium , Cambridge Universtiy Press, Cambridge (1998); 
T. Vicsek, Fractal growth phenomena, World Scientific 
(1989). 

K. De'Bell, A.B. Maclsaac, and J.P. Whitehead, Rev. 
Mod. Phys. 72, 225 (2000). 

P. Meakin and M. Muthukumar, J. Chem. Phys. 91, 3212 
(1989). 

N. Vandewalle and M. Ausloos, Phys. Rev E 51, 597 
(1995). 
R. Pastor-Satorras and J.M. Rubi, Phys. Rev. E 51, 5994 



(1995); Prog. Colloid. Polym. Sci. 110, 29 (1998); J. 

Magn. Magn. Mater. 221, 124 (2000). 
[7] P.M. Mors, R. Botet, and R. Jullien, J. Phys. A 20, L975 

(1987). 
[8] G. Helgesen, A.T. Skjeltorp, P.M. Mors, R. Botet, and 

R. Jullien, Phys. Rev. Lett. 61, 1736 (1988). 
[9] F. de los Santos, M. Tasinkevych, J.M. Tavares, and 

P.LC. Teixeira, J. Phys.: Condens. Matter 15, S1291 

(2003). 
[10] The relaxation of the dipole does not appear to be an 

essential ingredient of the algorithm, since the same D 

has been found in with and without dipole moment 

relaxation. 
[11] J.M. Tavares, M. Tasinkevych, F. de los Santos, and 

M.M. Telo da Gama, Mol. Phys. 101 1659 (2003). 
[12] P. Meakin, Phys. Rev. B 30, 4207 (1984); P. Meakin, J. 

Kertesz, and T. Vicsek, J. Phys. A 21, 1271 (1988); Z. 

Racz and T. Vicsek, Phys. Rev. Lett. 51, 2382 (1983). 
[13] P. Meakin and F. Family, Phys. Rev. A 34, 2558 (1986). 
[14] J.M. Tavares et al., unpublished. 



